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Abstract 

SLR technique has being used for determination of ERP during over twenty years. 
Most of results contributed to lERS are based on analysis of observations of Lageos 
1&2 satellites collected at the global tracking network of about 40 stations. Now 
5 analysis centers submit operative (with 2-15 days delay) solutions and about 10 
analysis centers yearly contribute final (up to 23 years) ERP series. Some statistics 
related to SLR observations and analysis is presented and analyzed. Possible problems 
in SLR observations and analysis and ways of its solution are discussed. 

1 Introduction 

Laser ranging to Earth artificial satellites (SLR) was initiated in 1964 after launch of the 
first geodetic/geodynamical satellite Beacon-B. 

Since that time satellite laser ranging (SLR) technique have being widely used for geody- 
namical and geophysical researches. The primary fields of investigations used SLR observa- 
tions are Earth rotation, maintenance of the Terrestrial Reference Frame, tectonic motion. 
Earth crust deformations, geopotential with its spatial and temporal variations, tides, move- 
ment of geocenter, support of satellite geophysical missions (such as satellite altimetry), 
global time transfer, and others. Detailed analysis of SLR contribution to Earth sciences 
can be found in (Tapley et al. 1993). 

After the launching of the Lageos satellite in May 1976, SLR became one of the main 
techniques for investigations of the Earth rotation, and during over twenty years SLR tech- 
nique have being used for determination of ERP. Most of results contributed to lERS are 
based on analysis of observations of Lageos 1&2 satellites collected at the global tracking 
network of about 40 stations. SLR provides high precision series of Xp, Yp, and LOD. Some 
analysis centers compute also UT that allow to densify Universal Time series in combination 
with VLBI data. 

Importance of ILRS as one of the main method to study the Earth led, naturally, to the 
establishment of the International Laser Ranging Service (ILRS) in 1998 that coordinates 
now scientific activity SLR (and LLR), chiefly in the framework of lAG and lERS projects. 

In this paper we will focus only on SLR contribution to investigation of polar motion in 
accordance with topics of the conference. 



*In: ASP Conf. Ser., v. 208, Polar Motion: Historical and Scientific Problems, Proc. lAU Coll. 178, 
Cagliari, Italy, Sep 27-30, 1999, Eds. S. Dick, D. McCarthy, B. Luzum, 2000, 267-276. 
^Current affiliation; Pulkovo Observatory, St. Petersburg, Russia 
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Figure 1: Contribution of SLR AGs to lERS yearly solutions. 

2 Contribution of SLR PM series to the lERS 

SLR observations of Lageos and ERP derived from these observations are available from 
May 1976. However, only in 1980-81 after significant improvement of range technique SLR 
ERP series achieved accuracy required for investigation of Polar motion (PM). Figure [Hsum- 
marizes information from lERS Annual Reports for 1978-1997 concerning use of submitted 
SLR series in lERS yearly solutions. One can see that beginning from 1983 SLR ERP series 
is regularly used for lERS combination. It's also seen that the Center for Space Research, 
University of Texas at Austin (GSR) provides most long-time spanned and stable PM series. 

At present 5 analysis centers submit operational solutions with 2-15 days delay, and 
about 10 analysis centers yearly contribute final (up to 23 years) ERP series. Most long 
series available for analysis are listed in Table [1] 

3 Accuracy of SLR PM series 

During long history of using SLR technique to study PM drastic improvement in range 
precision was achieved. It's interesting to see how this technology development affected 
accuracy of PM series obtained from SLR observations. Figure [2] shows differences between 
SLR series obtained at the GSR and the Institute of Applied Astronomy, St. Petersburg 
(lAA) and lERS combined solution EOP(IERS)G04. 

One can clearly see in Figure [2]that replacement of the first generation ranging equipment 
by the second generation one about 1980 and its further replacement by the third generation 
units leaded to significant improvement of accuracy of PM solution. However, in spite of 
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Table 1: Longest series. 



CSR 

GSFC 

DUT 

GAOUA 

lAA 

CGS 

BKG 



1976 
1980 
1983 
1983 
1983 
1985 
1988 



Table 2: Agreement of series 
with C04, mas. 



Solution Xr, 



GPS 
SLR 
VLBI 



0.11 0.11 
0.11 0.12 
0.16 0.15 



Table 3: Statistics of operational solu- 
tions based on USNO data for Sep 1998 
- Sep 1999. 
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Figure 2: Differences between SLR and lERS PM series. 
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improvement of ranging technique and implementation the fourth generation equipment, 
precision of PM series derived from SLR observations remains the same during last 10-12 
years. We will try to discuss the problem below. 

In 1998 lERS Central Bureau had derived final lERS combination in two step: on the 
first stage combined VLBI, SLR, and GPS series was computed and then these was used 
for final lERS combination. Table [2] copied from the 1997 lERS Annual Report shows that 
precision of SLR and GPS combinations is practically the same. Unfortunately, this is not 
so for operational solutions. 

Table [3] compiled on basis of gpspol.asc files produced by USNO along with Bulletin 
A issues contains statistics of series used for lERS Rapid Service. Of course, errors in 
PM obtained there depend on weighting applied to combined series, but it already clear 
that accuracy of individual SLR series is worse than VLBI and GPS ones, and delay of SLR 
contribution is much worse than GPS one. Nevertheless, one can see that at least two centers 
can produce operational solution with delay about 2 days. If all SLR analysis centers would 
provide such a solution on regular basis, it will allow to have, in principle, combined SLR 
solution with accuracy 0.1-0.15 mas and delay 1.5-2 days. 

4 Factors limited quality of the SLR ERP solutions 

During long time SLR was the one of the main methods of determination of PM and den- 
sification of PM series based on VLBI observations. However, lately SLR ERP series are 
inferior to GPS ones in quality (accuracy, density, delay of operational solution) of results. It 
would be very important to understand existing problems in SLR observations and analysis 
procedures and discuss possible ways to solve them. 

4.1 Observations 

The first problem with SLR observations is that SLR technique is an 'one-object' one. This 
means that, unlike GPS, SLR station can observe only one object in time. Hence, planning 
of observations and priority politics plays substantial role in acquisition of observations 
of satellites on which PM determination is based, especially keeping in mind that SLR 
observations are rather expensive and number of ranging units are limited. 

At present four operational satellites seem most suitable for investigation of Earth ro- 
tation (as well as tectonic motion and long-term temporal variations in various geophysical 
parameters) — two Lageos satellites and two Etalon satellites. Both Lageos and Etalon 
satellite was launched to long-time stable orbits and have a low area to mass ratio. Their 
description is presented in Figure [31 

Figure H] shows number of observations of these satellites during last 12 years and ILRS 
priorities. One can see that number of observations of Lageos satellites remains approxi- 
mately the same during these years and number of observations of Etalon satellite is too 
small to contribute seriously to analysis. 

Let us see in more details distribution of SLR observations in stations and time. Figure [5] 
shows distribution of observations in stations for the period from Sep 1983 till Aug 1999. 
Table H] contains list of stations contributed more than 2% of total number of observations 
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Lageos 



Etalon 





Satellite Satellite parameters 

D,in M,kg A/M,in'^2/kg 
Lageos 1,2 0.6 407 0.00070 
Etalon-1,2 1.3 1344 0.00098 



Orbital parameters 
a, km h, km T, h i, deg 
12300 5900 3.76 110,53 
25500 19100 11.26 65 



Figure 3: Satellites Lageos and Etalon. 



for whole period and during the last year (in parenthesis). 

The first impression is that stations are distributed quite uniformly along of satellite orbit. 
However, if we recall that SLR observations are weather-dependent, it would be evident that 
Australia - West Asia region need at least one more active station. Fortunately, several new 
stations was put into operation in China and Japan during the last years but they are still not 
active enough. So, problem of more uniform distribution of observations in stations evidently 
still exist. Besides, as usually, most of stations are located in the Northern hemisphere and 
at least two active stations would be very desirable in the Southern hemisphere (in South 
America and South Africa). 

Distribution of observations in months is more or less uniform, with small decrease in 
November and December (sharp dip in Christmas days!). But distribution in days of week 
(Figure [6]) is not so good. One can see that number of observations on several most active 
stations is about twice lesser than during workdays. This is statistics for whole interval from 
Sep 1983 till Aug 1999. For the last year distribution is a little bit more uniform - maximum 
number of observations was made in Wednesday (16.2%), minimum in Sunday (11.6%), but 
still far from ideal. 



4.2 Analysis 

Though, in my opinion, the only real way to improve quality of SLR ERP series is to increase 
number of observations and, especially, number of satellites involved in determination of 
ERP, analysis strategy can also be advanced. 

Modern tendency is to use combined VLBI, GPS, SLR series for final lERS combination. 
Unfortunately, PM series produced by SLR analysis centers are not unified enough which 
makes difficult their comparison and combination. The main problems in individual SLR 
series are: 
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Figure 4: Number of observations (normal points), thousands (at the left) and ILRS priorities 
(at the right) of satellites. 



• Sparse time series (3-5 days for most solutions, 1 day needed). 

• Large delay of operational solutions (5-10 days for most solutions, 1-2 days needed). 

• Non-unified TRS. 

• Dependence on a priory values of EOP. 

Analysis of methods used for determination of ERP from SLR observations shows that not 
all possibihties for improvement and unification of SLR ERP series are exhausted. What 
can we do to make SLR series more efficient? Existing experience shows that significant 
improvement could be achieved if we would: 

• Produce daily EOP series. 

• Produce operative solution with delay about 2 days. 

• Advance strategy for determination UT. 

• Use the same TRS for all solutions (ITRF seems most suitable for this purpose) . 



6 




Figure 5: Distribution of observations in stations. 



Table 4: Distribution of observations in stations. 
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Total about 130 stations including mobile occupations 

5 stations provided 38% of total number of observations 
7 stations provided 48% of total number of observations 
10 stations provided 62% of total number of observations 
15 stations provided 75% of total number of observations 
22 stations provided 85% of total number of observations 
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Figure 6: Distribution of observations in days of week for all stations (left) and selected 
stations (right). 



An analysis of these problems and some proposals on their solution was presented in 
(Malkin 1998). Below proposed methods tested in the lAA during last three years are 
briefly described. 

The first serious problem is how to compute daily ERP series without loss of accuracy. 
There are two main ways to compute daily series. The first is to use daily arc solutions (i.e. 
use one day arcs). But this leads to degradation of accuracy of solution. The second one is to 
use overlapping arcs with appropriate length and one day shift between arcs. Such a solution 
produces dependent results for neighbor arcs, but significant improvement in precision can 
be considered as sufficient compensation. Table |5] shows accuracy of ERP solution with 
various length of arc. 

So, experience shows that using overlapping arcs yields more precise ERP series than 
using one day arcs. 

It would be very desirable to have strictly daily series for all analysis centers, i.e. series 
with epochs equal to O'^. This can be achieved by two methods - determination of Pole 
coordinates with their rates and a posteriori interpolation to O'^. Numerous test showed that 
including Pole coordinates rates into parameter model lead to small degradation of accuracy 
of PM, evidently due to degradation of normal system matrix condition. A posteriori inter- 
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Table 5: Dependence of errors in ERP on length of arc (RMS of EOP(IAA)L ~ 
EOP(IERS)C04 after removing trend). 

Length of arc, days 



X, 0.001" 0.40 0.23 0.20 0.18 0.16 

y, 0.001" 0.34 0.22 0.21 0.18 0.15 

LCD, 0.0001s 0.44 0.16 0.15 0.12 0.09 



Table 6: Effect of interpolation of SLR ERP series (RMS of EOP(IAA)L - EOP(IERS)C04 
after removing trend). 

EOP Series 

Raw Interpolated 

X, 0.001" 0.154 0.155 

y, 0.001" 0.171 0.171 

LCD, 0.0001s 0.154 0.168 



polation of ERP to epochs gives better result. Linear interpolation was found as most 
suitable for this purpose. 

Table |6] shows dependence of results on effect of interpolation of raw ERP series. One 
can see that use of hnear interpolation don't lead to visible degradation of accuracy. Special 
tests showed that strictly daily series obtained with interpolation is more accurate that series 
obtained with including Pole coordinates rates in parameterization. 

It should be mentioned also that producing daily solution (independently which method 
id used) provide operational solution with steady delay about 2 days, which solves the second 
problem mentioned above. 

Other serious problem is determination of UT from SLR (and other satellite observations). 
It is well known that UTl cannot be separated from longitude of node of satellite orbit during 
parameter solution. Three methods are being used to solve this problem: 

• Fixing longitude of node during parameter solution (usually, during last iteration). 

• Analysis of node longitude series, forecasting it and use predicted values for operational 
solution. 

• Integrating LOD series to obtain independent free-running UT series with its possible 
correction for high-frequency variations from comparison with VLBI series. 

Evidently, only the latter method can provide (in principle) independent result. Since 
that is not a subject of this paper, we will not stay on detailed analysis of this problem. 
However, It is worth to mention here that significant improvement of SLR UT series is also 
impossible without increasing of number of satellites involved in determination of ERP. 
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At last, use the same terrestrial reference frame for all SLR solutions seems evident to 
achieve uniform solutions for combination. Use of ITRF as terrestrial reference frame for 
by all analysis centers for their SLR solutions provides more homogeneous series for SLR 
combined solution. 

Realization of this or alternative analysis strategy could provide more uniform, accurate 
and operative SLR ERP series. After that, combining of all submitted series to final ILRS 
SLR product seems reasonable and useful for further use for lERS and other purposes. 

We have not mention here such a serious problem as dependence of SLR ERP results 
on a priori values. This is worth to perform special investigation for each method used for 
computation of ERP at various analysis centers. 

5 Conclusion 

Satellite laser ranging technique made and make a very valuable contribution to Earth dy- 
namics. In particular, very valuable contribution was made in investigation of PM. During 
many years SLR was one of the main methods of determination of polar motion and main 
method of densification of ERP series obtained with VLBI. At present 5 analysis centers 
submit operational solutions with 2-15 days delay, and about 10 analysis centers yearly 
contribute final (up to 23 years) ERP series. 

However, due to principal peculiarity of this method (relatively expensive experiment, 
limited number of units, lack of capability of multi-satellite ranging, etc.) quality of SLR 
ERP data remains the same during the last decade in spite of ranging precision improved by 
a factor of a thousand from a few meters to a few millimeters since the first SLR experiments 
(and by factor of about 10 during last 10-12 years). This leads to decreasing of weight of 
SLR solutions in the combined lERS EOP series. 

It is evident that only substantial increasing of number of observations and satellites 
involved in investigation of Earth rotation and improvement of distribution of observations 
in stations and time can help in improvement of the SLR ERP accuracy. However, capacity 
of tracking network is practically exhausted. 

On the other hand, in spite of GPS provides determination of ERP with impressive accu- 
racy and delay, SLR results are very important for combined lERS solution for improvement 
of systematic accuracy of the final lERS product. Analysis of precision of individual SLR 
series shows that its combination can provide combined SLR much more accurate and rapid 
series. To achieve highest accuracy of combined SLR ERP product is necessary to solve 
problems discussed in section 14.21 

Realization of this opportunities by ILRS would be very important for investigation 
of Earth rotation because allow to save independent method of determination of PM and 
velocity of the Earth rotation. 
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